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ABSTRACT 

I propose an explanation to the finding that star formation and visible fila- 
ments strong in Ha emission in cooling flow clusters occur only if the minimum 
specific entropy and the radiative cooling time of the intracluster medium (ICM), 
are below a specific threshold. The explanation is based on the cold feedback 
mechanism. In this mechanism the mass accreted by the central black hole orig- 
inates in non-linear over-dense blobs of gas residing in an extended region of the 
cooling flow region. I use the criterion that the feedback cycle period must be 
longer than the radiative cooling time of dense blobs for large quantities of gas 
to cool to low temperature. The falling time of the dense blobs is parameterized 
by the ratio of the infall velocity to the sound speed. Another parameter is the 
ratio of the blobs' density to that of the surrounding ICM. By taking the values 
of the parameters as in previous papers on the cold feedback model, I derive an 
expression that gives the right value of the entropy threshold. Future studies will 
have to examine in more detail the role of these parameters, and to show that the 
observed sharp change in the behavior of clusters across the entropy, or radiative 
cooling time, threshold can be reproduced by the model. 

Subject headings: galaxies: clusters: general — cooling flows — galaxies: active 
— intergalactic medium — X-rays: galaxies: clusters 



1. INTRODUCTION 



In the last decade it has become clear that the intra-cluster medium (ICM) in cooling 
flow (CF) clusters of galaxies and CF galaxies must be heated, and the heating process 
should be stabilized by a feedback mechanism (see review by McNamara & Nulsen 2007). 
However, in many cases the heating cannot completely offset cooling (e.g., Wise et al. 2004; 
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McNamara et al. 2004; Clarke et al. 2004; Hicks & Mushotzky 2005; Bregman et al. 2006; 
Salome et al. 2008) and some gas cools to low temperatures and flows inward (e.g., Peterson 
& Fabian 2006). The mass inflow rate is much below the one that would occur without 
heating, and the flow is termed a moderate cooling flow (Soker et al. 2001; Soker & David 
2003; Soker 2004). 

One of the likely ingredients of the moderate CF model is a cold feedback process 
(Pizzolato & Soker 2005; Soker 2006; Pizzolato 2007). In the cold-feedback model the 
mass accreted by the central black hole originates in non-linear over-dense blobs of gas 
residing in an extended region of r ~ 5 — 30 kpc; these blobs are originally hot, but then cool 
faster than their environment and sink toward the center. The mass accretion rate by the 
central black hole is determined by the cooling time of the ICM, the entropy profile, and the 
presence of inhomogeneities (Soker 2006). Most important, the ICM entropy profile must be 
shallow for the blobs to reach the center as cold blobs. This accretion process is different from 
the commonly assumed accretion mode in feedback models where the black hole accretes hot 
gas (hot feedback mechanism) from its vicinity, e.g., via a Bondi-type accretion flow (e.g. 
Churazov et al. 2002; Nulsen 2004; Omma & Binney 2004; Chandran 2005; Croton et al. 
2006; Balmaverde et al. 2008). 

Among other things, the cooling of the ICM to low temperatures is inferred from star 
formation (e.g., Rafferty et al. 2008), and from cool filaments via their Ha emission (Ha 
filaments). In a recent paper Cavagnolo et al (2008) found that almost all clusters with 
strong H-alpha emission have a central entropy of K = kTrie 2 ^ 3 < 30 keV cm 2 ; here k is the 
Boltzmann constant, T is the ICM temperature, and n e its electron density. This relation 
might be better presented as a relation between star formation (SF) and the radiative cooling 
time: no star formation is seen if the radiative cooling time is r coo i > 5 x 10 8 yr or the entropy 
is K > 30 keV cm 2 (Rafferty et al. 2008). 

Voit et al. (2008) suggest that this K— Ha (or SF— r coo i relation) results from a compe- 
tition between radiative cooling of cooler regions and heat conduction into these regions. In 
high entropy ICM the heat conduction manages to prevent thermal instabilities that poten- 
tially could form Ha filaments. For the cooling function dependance on temperature Voit 
et al. (2008) assumed A oc T 1 / 2 . How ever, for most clusters the temperature in the inner 
region is T < 3 keV, and for a composition above half a solar the cooling function is basically 
constant in the range ~ 1.7 — 3 keV. This introduces a very small change in the numerical 
value, and is not significant. But it does show that the dependance on some power of the en- 
tropy is an approximate one. The same holds for the explanation proposed here (see section 
3). More significant is their assumption that the size of the cooling region is of the order 
of the radius of the cluster at the location of the region. The instabilities at such scale are 
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likely to be very fragmented (fractal), such that the heat conduction is more efficient than 
what they assume (because the temperature gradient is larger across the smaller dimension 
of a filament, and the surface area for heat conduction is much larger). Also, if in some 
regions the magnetic field lines are closed on themselves on a much smaller scale than the 
distance to the cluster center r, then the heat conduction is limited to within magnetic flux 
tubes (Soker et al. 2004; Soker 2004). Another point of concern is discussed in section 3. 

A detail study of the competition between radiative cooling of cool filaments, at a 
temperature of ~ 10 4 K, and heat conduction from the ICM is presented by Nipoti & Binney 
(2004). They find that cool filaments survive only if the ICM density is high and the ICM 
temperature is low, i.e., has a low entropy. However, Nipoti & Binney (2004) argue that the 
filaments do not originate from the hot ICM as Voit et al. (2008) consider, but rather the 
cool gas originates from other sources, e.g., AGN activity and external infall. 

Voit et al. (2008) note that the radiative cooling rate of the ICM depends on entropy, 
and that the K— Ha relation might be alternatively related to the short radiative cooling 
time of the low-entropy ICM, hence better termed SF— r cool relation (Rafferty et al. 2008). 
In this paper I use the competition between radiative cooling time and the infall time of 
dense blobs as developed for the cold feedback mechanism by Pizzolato & Soker (2005), to 
explain the K— Ha (SF— r coo i) relation. 



In the cold feedback mechanism cold dense blobs fall toward the center of the cluster 
and feed the central AGN. Say a region in the cluster is perturbed, and blobs start to cool. 
As the first blobs fall and reach the center, an AGN outburst will heat the perturbed region, 
and prevent further cooling. The condition for large quantities of gas to cool is for the 
feedback cycle to be longer than the radiative cooling time 



The cooling time of a blob with a density of 5 times the ambient density, hence a 
temperature 5 times lower, is given by 



2. 



THE COLD FEEDBACK EXPLANATION 




(1) 




(2) 



where the total, electron, and proton densities of the ICM, respectively, and 

T is the ICM temperature. In equation (j2j) the increase in the cooling function A and in 
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the density with decreasing temperature of the cooling blob was accounted for. For the 
relevant temperature range here, 1.5 keV < kT < 4 keV, and for solar composition the 
cooling function is A ~ 2 x 10~ 23 erg cm 2 s _1 (Gaetz et al. 1988). 

For the feedback cycle time I assume the following. The dense blob falls with its terminal 
velocity v t . Pizzolato & Soker (2005) studied the fall of dense blobs through the ICM, and 
found the terminal velocity to be 



100 pc/ vio- 8 cms- 2 / V 3 

where a is the radius of a spherical falling blob, g is the cluster's gravitational field, and 
the same scaling as theirs was used here. In this paper I scale the blob size and density 
contrast with the typical values used by Pizzolato & Soker (2005). I incorporate the radius 
of the blob and 5 into a parameter rj, such that the average inflow velocity of the blob is 
7]C S = i](5kT ' l?> iimn) 1 ! 2 \ where C s is the sound speed of the ICM. For a cluster temperature 
of 2 keV, a velocity of 50 km s _1 corresponds to rj = 0.07. The feedback cycle period is then 

r 

tcyc = ^r- (4) 

Inserting equations (E])-(j4j) in condition (TjQ) gives the condition for the formation of a 
large quantities of cold gas 

r > I) e = r,C " kT 



5 2 n e n p A 



37 (2 x 10- 23 erg cm 2 s" 1 J (lO keV cm 2 J kpC ' ^ 

where K = kTrie 2 ^ is the entropy. The above condition implies that for blobs to cool within 
the cooling radius, the entropy must be lower than some threshold. 

The power of the entropy in equation ([5]) is an approximate one. The value of 3/2 was 
obtained under the assumption that the infall speed of dense blobs is proportional to the 
sound speed, and on that the cooling gas is in the temperature range where the cooling 
function does not depend on temperature. The last assumption is a good approximation 
for most clusters near the threshold K ~ 30 keV cm 2 . In any case, in the cold feedback 
mechanism what matters is the radiative cooling time. This is also seems to be the case 
from observations (Rafferty et al. 2008). 

Condition ([5]) is sensitive to the value of 5 and rj. Here I simply took the same value used 
by Pizzolato & Soker (2005). In obtaining rj I assumed a blob size of ~ 100 pc, as the scaling 
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used by them, and for 5 I took the value of an unstable blob used by them. The parameter 
r) has in it the dependance on the blob radius as a 1//2 , and on the density contasrt as 5 1 / 2 . 
Thereofre, condition (jSJ) practically depends on a 1 / 2 and S 3 ^ 2 . The limit value of 1.1 kpc in 
condition (J3J) is similar to the one obtained by Voit et al. (2008). With their typical value 
of the heat conduction suppression factor f c = 0.2 their equation (2) has a limit of 1.8 kpc. 
However, they use a lower value for the cooling function (because they take the dependance 
to be A oc T 1 / 2 ). If I take a value for A as they take, then the two coefficient are almost 
equal. Still, the sensitivity of condition §5§ to the parameters r\ and 5 is a somewhat weak 
point of the proposed explanation. 

The cold feedback mechanism can account also for the finding that the low entropy 
clusters have stronger radio emission (Donahue et al. 2005; Cavagnolo et al. 2008). The 
same source of blobs that lead to the formation of Ha filaments, will feed the AGN activity. 
Cavagnolo et al. (2008) noticed that below the threshold K < 30 keV cm 2 there is no 
correlation between radio power and the central entropy. They tentatively speculated that 
this lack of correlation hints that cold-mode accretion (Pizzolato & Soker 2005; Hardcastle 
et al. 2007) might be the dominant process feeding the AGN. 

3. DISCUSSION AND SUMMARY 

In this paper I showed that the findings of Rafferty et al. (2008) that no star formation 
(SF) is seen if the radiative cooling time is r coo i > 5 x 10 s yr, and of Cavagnolo et al. 
(2008) that almost all clusters with strong H-alpha emission have a central entropy of K = 
kTn~ 2/3 < 30 keV cm 2 , might be explained by comparing the cooling time of dense blobs 
in the ICM (eq. H]) with the response time (cycle time) of the AGN feedback heating. This 
time is taken to be equal to the fall time of the dense blobs to the center (eq. H]). Many 
dense blobs will cool to low temperatures if the response time is longer than their cooling 
time (eq. [p. This leads to equation ([5]) which is the main result of this paper. It shows that 
for blobs to cool within the cooling flow (CF) radius, the entropy must be lower than some 
threshold. 

The feedback heating of the ICM in this model is maintained by dense blobs that are 
accreted by the central black hole and originate in non-linear over-dense blobs of gas residing 
in an extended region of r ~ 5 — 30 kpc; this is the cold feedback mechanism (Pizzolato & 
Soker 2005; Soker 2006). 

Criterion ([5]) can be compared with the one derived by Voit et al. (2008) 
r > 4(A7l0 keV cm 2 ) 3//2 / c 1//2 , where f c is the suppression factor of the heat conductivity, 
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taken to be f c ~ 0.1 — 1 in their model. The similarity of the expressions are interesting. They 
have the same dependance on the entropy a similar numerical factor, and some parameters. 
Voit et al. (2008) consider their expression to satisfactorily explain their observations. The 
similarities between the two conditions show that the cold feedback explanation proposed 
here should be considered in future studies as well. In particular, the result is satisfactory 
if we consider the crude derivation performed here, and the inhomogeneous nature of the 
ICM, where large regions with lower entropy than the average exist. The dense blobs will 
be more likely to develop in these regions. 

The value of f c ~ 0.1 — 1 and the large thermally unstable regions considered by Voit 
et al. (2008) implies that heat conduction is globally important in the CF region. However, 
models based only on global heat conduction are unstable and require fine tuning (Bregman 
& David 1988; Soker 2003; Kim & Narayan 2003). The explanation of Voit et al. (2008) 
requires a stabilizing mechanism. In a recent paper Guo et al. (2008) built a feedback model 
based on both heat conduction and AGN heating. They showed that the AGN heating 
stabilizes the feedback mechanism. 

There is a price payed to achieve the stability in the model of Guo et al. (2008). (1) The 
inner boundary of their simulation is at r = 1 kpc. They require that all the mass that enters 
the r = 1 kpc sphere will reach the central BH, and the mechanical energy of the launched 
jet be L agn = eM(l kpc)c 2 , with e ~ 0.1 — 0.3. This is an extremely efficient conversion of 
mass inflowing at r = 1 kpc to mechanical energy of the jets. (2) In their stable model of 
A1795 (their model A3) the entropy decreases with radius, such that the model is unstable 
to convection. This might not be a big concern, as by using a different AGN heating 
prescription the negative entropy might be removed (Guo, F., private communication 2008). 
(3) In their best models of A2199 (model B3) and of A2052 (model C3) the response time 
of the AGN is too long. Using their values of the mass inflow rate and density at r = 1 kpc, 
the inflow velocity Vf(l kpc) can be calculated. The respond of the accreting black hole to 
changes in the ICM is on a time scale of r r = 1 kpc/f/. I find r r = 7 Gyr for model B3 and 
r r = 8 Gyr for model C3. In both cases the response time r r is about an order of magnitude 
longer than the cooling time in the center of these clusters (0.6 and 1.1 Gyr, respectively). 
This shows that the AGN heating has no time to respond to changes in the thermal state of 
the ICM. It is possible that this problem will disappear if an episodic AGN heating will be 
applied (Guo, F., private communication 2008). 

The explanation proposed here for the K— Ha relation (or SF— r coo i relation) has to 
overcome some difficulties as well. In the first step it will be required to show that reasonable 
values for the two parameters used in the model (<5 and rf) can lead to the sharp transition 
from bright Ha clusters, for K < 30 keV cm 2 , to clusters with no, or very weak , detection 
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of Ha emission. Here I note that the values of 77 and S used here are the typical values used 
by Pizzolato & Soker (2005). If I insert the threshold value of K = 30 keV cm 2 found by 
Cavagnolo et al. (2008), and 5 ~ 2.5 — 3 used by Pizzolato & Soker (2005), in equation fl5j) 
I find the numerical value there to be ~ 13 — 6 kpc. This is inside the range of r ~ 5 — 30 kpc 
considered by Pizzolato & Soker (2005) to be the region where dense blobs originate. In a 
later step, numerical simulations of dense blobs in the ICM should show that the formation 
of cold clouds that are the source of Ha emission depends on entropy in the observed way. 

The main difference between the explanation proposed here and that of Voit et al. (2008) 
is the importance of the global heat conduction. In recent years the debate on the question 
of whether global heat conduction in CF clusters is important or not has intensified. If that 
debate is resolved toward a significant global heat conduction, then the explanation proposed 
here must be ruled out. If it will turned out that global heat conduction is not important, 
then the explanation of Voit et al. (2008) must be ruled out. In the cold feedback mechanism 
cold blobs feed the central BH. Therefore, in CF clusters with strong AGN activity, some 
cold blobs exist close to the center. Therefore, the prediction here is that very low level of 
Ha emission might be detected even in CF clusters with high entropy. 

I thank Mark Voit, Brian McNamara, Fabio Pizzolato, and Fulai Guo for helpful com- 
ments. This research was supported in part by the Asher Fund for Space Research at the 
Technion. 
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